Exotic Mott phases of the extended t-J model on the checkerboard lattice at 

commensurate densities 
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Coulomb repulsion between electrons moving on a frustrated lattice can give rise, at simple com- 
mensurate electronic densities, to exotic insulating phases of matter. Such a phenomenon is il- 
lustrated using an extended t-J model on a planar pyrochlore lattice for which the work on the 
quarter- filled case cond-mat/0702367 is complemented and extended to 1/8- and 3/8-fillings. The 
location of the metal-insulator transition as a function of the Coulomb repulsion is shown to depend 
strongly on the sign of the hopping. Quite generally, the metal-insulator transition is characterized 
by lattice symmetry breaking but the nature of the insulating Mott state is more complex than 
a simple Charge Density Wave. Indeed, in the limit of large Coulomb repulsion, the physics can 
be described in the framework of (extended) quantum fully-packed loop or dimer models carrying 
extra spin degrees of freedom. Various diagonal and off-diagonal plaquette correlation functions are 
computed and the low-energy spectra are analyzed in details in order to characterize the nature of 
the insulating phases. We provide evidence that, as for an electronic density of n=l/2 (quarter- 
filling), the system at n = 1/4 or n = 3/4 exhibits also plaquette order by forming a (lattice 
rotationally-invariant) Resonating-Singlet-Pair Crystal, although with a quadrupling of the lattice 
unit cell (instead of a doubling for n = 1/2) and a 4-fold degenerate ground state. Interestingly, 
qualitative differences with the bosonic analog (e.g. known to exhibit columnar order at n=l/4) 
emphasize the important role of the spin degrees of freedom in e.g. stabilizing plaquette phases 
w.r.t. rotational symmetry-breaking phases. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg, 74.20. Mn, 71.10.Fd 



A. Introduction 

Frustrating antiferromagnetic (AF) interaction in 
quantum magnets^ can give rise to various exotic spin 
gapped quantum disordered phases which could be real- 
ized in a number of fascinating materials such as Kagome, 
spinels or pyrochlore materials. Among these phases, 
a Valence Bond Crystal (VBC), the so-called "plaque- 
tte phase" , which spontaneously breaks lattice symmetry 
(while preserving rotation symmetry), is of particular in- 
terest and has been identified^ in e.g. the Heisenberg 
model on the checkerboard lattice, a two-dimensional 
(2D) lattice of corner sharing tetrahedra (i.e. the 2D ana- 
log of the three-dimensional pyrochlore lattice as shown 
inFig.HD. 

Sofar, itinerant correlated fermions on frustrated lat- 
tices have been poorly explored and most investigations 
deal with weakly doped Mott insulators^. A bit more 
is known for bosonic systems (which can be mapped on 
quantum spins under magnetic field) for which a larger 
set of numerical techniques is available (like Quantum 
Monte Carlo, QMC). On the triangular lattice, hard-core 
bosons with nearest-neighbor (NN) repulsive interaction 
exhibit a rich phase diagram^ with charge ordering at 
commensurate 1/3 or 2/3 fillings as well as a supersolid 
(i.e. a phase with both superfluid and charge orders) 
under light doping. Recently, hard-core bosons with NN 
repulsion have also been investigated on the checkerboard 
lattice at 1 /4- filling where a weakly first-order superfluid- 
insulator transition has been detected when increasing 



Coulomb repulsion^. Interestingly, in the large repul- 
sion insulating phase, evidences have been found in favor 
of lattice symmetry breaking associated to plaquette or- 
dering similar to the one of the Quantum Dimer Model 
(QDM) on the square lattice-. A similar plaquette phase 
has also been identified at half-filling in the same model 
and on the same lattice using an effective Hamiltonian 
valid at large repulsion^,. Such a fascinating behavior is in 
fact due to the partial particle localization occurring via 
"ice-rule" -like constraints enforced by the minimization 
of the global Coulomb energy, a "classical" frustrated 
minimization problem. It is crucial here that such con- 
straints preserve a macroscopic degeneracy of the clas- 
sical configuration manifold. Indeed, for quarter-filling 
the Hilbert space is given by the ensemble of fully packed 
dimer coverings (as for the QDM, see Ref. 0) and a fully 
packed loop representation holds at half-filling^. Quan- 
tum fluctuations which occur within these subspaces via 
second-order kinetic processes can subsequently lead to 
exotic symmetry-breaking states. Here the basic ingredi- 
ent is the interplay between the frustrated nature of the 
lattice and the Coulomb repulsion (the Ising interaction 
in the "spin language"). 

The investigation of the fermionic analog, i.e. fermions 
moving on frustrated lattices and subject to (short range) 
Coulomb repulsion at commensurate densities, is more 
involved due the famous "fermion sign problem" that 
prohibits large scale QMC simulations. However, lim- 
ited progress have been made recently thanks to the con- 
struction of effective models describing the low-energy 



2 



physics at large Coulomb repulsion. These models (or 
simple extensions of them) might have known GS in some 
very special limits (analogous to the Rokhsar-Kivelson 
point of Ref. ^) and, otherwise, are tractable by weak 
coupling expansion around a periodic pattern of appro- 
priately chosen clusters- and by Exact Diagonalisation 
(ED) techniques. In the case of spinless fermions on 
the checkerboard lattice, the dynamics within the con- 
strained Hilbert space is introduced via 6-sites ring pro- 
cesses on hexagons involving 3 particles simultaneously 
and gives rise to a rich phase diagram— and, under some 
conditions, to particle fractionalization^^. The case of 
spinful fermions which has been recently investigated at 
and around quarter-fillingi^ involves a second-order ki- 
netic process on the "empty squares" (i.e. the void pla- 
quettes) of the planar pyrochlore lattice. Interestingly, 
it was shown that the quarter-filled GS is a new type 
of plaquette insulating phase, the Resonant Singlet Pair 
Crystal (RSPC): qualitatively, within a given sublattice 
of the empty squares (the empty squares of the checker- 
board lattice of Fig.[T]being decomposed into a "checker- 
board" pattern of two sublattices) , electron singlet pairs 
formed on the square diagonals resonate between the 
two possible configurations. It was shown that such a 
plaquette phase is remarkably robust w.r.t. perturba- 
tions (like diagonal empty square repulsion) contrary to 
the bosonic case discussed earlier for which a plaquette 
phase also exists but with a much more limited extension 
in parameter-space. In addition, it has been argued that 
lattice symmetry breaking would survive under light dop- 
ing in the hole-paired region, giving rise to a fermionic 
analog of a supersolid^. It should be noticed that, be- 
sides electronic systems on frustrated lattices with large 
Coulomb interactions, cold atoms in optical lattices offer 
also a fantastic playground for such ideas since models of 
the form discussed in this paper can in principle be real- 
ized by the appropriate tuning of atomic interactions^^. 

In this paper, we largely extend the investigation of 
correlated fermions on the checkerboard lattice initiated 
in Ref. [l^ by considering, besides ?i = 1/2 (quarter- 
filling), other electron fillings like n — 1/4 (1/8-filling) 
and n — 3/4 (3/8-filling) potentially interesting to stabi- 
lize new insulating phases. We study the metal-insulator 
transition within the extended t-J model for increas- 
ing Coulomb repulsion and we introduce then (follow- 
ing Ref. [12) an effective model describing the insulating 
state of matter. In order to characterize the various in- 
sulating phases, (i) we analyze the structure of the low- 
energy spectrum in the light of group theory predictions 
for symmetry-breaking phenomena and (ii) we compute, 
for the first time in the case of itinerant systems, diagonal 
and off-diagonal plaquette correlation functions. 

Our starting point is the (fermionic) extended 
Hubbard-like model including nearest neighbor repulsion 
V on the checkerboard lattice. Here, we consider the 
limit where the on-site repulsion U is very large (and 
the largest energy scale), thus forbidding double occu- 
pancy via the Gutzwiller projector Vg and leading to a 
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FIG. 1: (Color on-line) Left: the pyrochlore lattice is com- 
posed of "crossed plaquettes" and "empty squares" . The two- 
site unit cell (dashed square) and the X and Y axes considered 
in this study are shown. The C4v point group is defined w.r.t 
the center O of an empty square. Right: Path T ~ M — K — F 
in the first Brillouin zone considered in this study. 

t-J-V model description: 

n = -tJ2'PG{clc^„ + h.c}jVG + Vj2^^^J+Hj, 
Hj ^'J^{\-^i-^j\ n^n, , (1) 

where J is the AF exchange constant. Hereafter, when 
dealing with this model and if not specified, t = 1 sets 
the energy scale. Exact Diagonalisations (ED) have been 
performed on a periodic a/32x -\/32 cluster {N = 32 sites) 
using full translation, point group (see Fig. [T]) and time- 
reversal {Sf —Sf) symmetriesi^. 

B. Metal-Insulator transition in the t-J-V model 

In the following, we focus on the special electron den- 
sities, n = fc^, where the integer fc is 1, 2 or 3, for which 
we expect an insulator for sufficiently large ratio V/\t\. 
The metal-insulator (MI) transition (for increasing V) we 
are interested in is not directly connected to any nesting 
properties of the Fermi surface but it is rather induced by 
the Coulomb potential itself (hence it appears for both 
signs of t). To understand it qualitatively, let us first 
assume that the NN repulsion V is the largest energy 
scale in Hamiltonian ([T|). In that case, the Coulomb en- 
ergy is minimized by fulfilling the "ice rule" constraint 
of exactly k electrons in each tetrahedra. This condition 
still preserves an extensive (e.g. macroscopic) manifold of 
states^ii^iiiii^. Although quantum fluctuations occur be- 
tween such configurations (in second order in i), the GS 
is an insulator since no charge can be driven between two 
distant locations in the system without violating the local 
constraints. Because of the huge (i.e. macroscopic) de- 
generacy of the classical problem {V = oo), it is clear that 
the GS should be different from a simple CDW where par- 
ticles are localized at fixed locations, and should exhibit 
a more subtle and richer nature as seen later on. 

We have investigated numerically the MI transition as 
a function of V by computing the charge structure factor 
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plaquettes) occurs for arbitrary small interactionlS while 
a different insulating phase is stabilized at larger V as dis- 
cussed already in Ref. [H and in the following sections. 
An Insulator-to-Insulator transition is then expected in 
contrast to n = 1/4 and n = 3/4 (for both signs of t) 
and in contrast to n = 1/2 with t < 0. In all the latter 
cases, a single MI transition occurs as depicted above for 
n = 1/4, a density for which a numerical calculations was 
feasable. 



Effective model and mapping onto Quantum 
Dimer or Loop models 



FIG. 2: (Color on-line) Density structure factor in the GS of 
the t-J-V model for n = 1/4 and J — 0.2 calculated by ED on 
a 32-site cluster, (a) Structure factor at momentum q = (tt, 0) 
or (0, tt) vs T^/|t|; (b,c) Structure factors as a function of mo- 
mentum following the path in the BZ shown in Fig. [T] Open 
(dashed) symbols correspond to positive (negative) t. The 
points labelled as "V » 1" are obtained with the effective 
model with W = 0. 

^(q) = *5xp (iq ■ R)-/V(|R|), the Fourier transform 
of the (equal-time) density-density spatial correlation at 
distance i? = |R|, 

I ' 

where the sum (over N /2 terms, N being the number of 
sites) is restricted to sites i and j such that jr^ — | = i? 
and R has integer coordinates in terms of the unit cell 
lattice vectors of Fig. [T] (hence i and j belong to the same 
sublattice). Note that the disconnected part has been 
subtracted for convenience. As shown in Fig. [2] for n — 
1/4 and J/t = 0.2, characteristic "Bragg" peaks at q = 
(tt, 0) and (0, tt) appear for increasing V signaling some 
form of charge ordering characteristic of a MI transition. 
However, the complexity of the insulating phase will be 
revealed by the simultaneous emergence of other types of 
correlations, such as plaquette correlations to be discussed 
later on, suggesting a more subtle and exotic GS than a 
simple CDW. We also note that the occurrence of a MI 
transition is of course independent on the sign of the 
hopping t. However, our data show that the rise of the 
Bragg peak is more abrupt for negative t and we cannot 
exclude the MI transition to be located at y = in that 
case (which might be related to the existence of a flat 
band at the bottom of the tight-binding band) . However 
the finiteness of the cluster does not enable to address 
the issue of the order of these transitions. 

Similar behaviors are of course expected for n = 1/2 
and n = 3/4 although a direct calculation with Hamil- 
tonian ([T]) is beyond the power of available computers. 
Note however that for t > perfect nesting of the non- 
interaction Fermi surface occurs at n = 1/2; hence an 
instability towards a VBC (with non-equivalent crossed 



The ice rule constraint (in the insulator at large- V) 
leads to simple interesting connections to quantum fully- 
packed loop or dimer models^. Indeed, by associating 
to each electron a dimer joining the centers of the two 
corner-sharing tetrahedra it belongs to, a configuration 
atn = l/4orn = 3/4(n = 1/2) can be represented 
as a dimer (loop) conhguration. In contrast to bosonic 
analogs^, each dimer now carries a spin-1/2 degrees of 
freedom. To fully characterize the insulating phase, fol- 
lowing Ref. m, we consider its effective Hamiltonian act- 
ing within the constrained Hilbert space, Ti = Hq + Hj 
with : 

Ha = -i2^Pn(s), (3) 

s 

Pais) = (cItS'i ~^uS't) (cfei^iT -CfcT^u) (4) 

where t2 = ^ and the summation is over the void pla- 
quettes labelled by s of the checkerboard lattice. The 
sites of s are ordered as i,k,j and I in the clockwise (or 
anti-clockwise) direction. Pais) is a second-order pro- 
cess which preserves the ice rule. It acts on two electrons 
forming a singlet bond on one of the two diagonals of s, 
and rotates the bond by 90° degree. 

D. Description of the expected ordered phases 

The kinetic processes of Eq. ^ favor singlet electron 
pairs resonating on void plaquettes (empty squares) and 
could lead to long-range plaquette-plaquette correlations. 
Fig. [31[a) shows pictorially the expected Resonant Singlet 
Pair Crystals exhibiting such a plaquette ordering. For 
n = 1/2 it was indeed shown that the q = (tt, tt) RSPC 
is stable- as long as J/t2 < 1.5. Note that the GS 
of the RSPC is two-fold degenerate at n = 1/2 while 
one expects a four-fold degenerate GS at n — 1/4 and 
n = 3/4 in connection with q — (tTjTt), (tTjO) and (0,7r) 
plaquette order as seen later. In addition, the average 
charge density of the n — 1/2 RSPC is uniform while 
q = (tt, 0) and (0, tt) (charge density) Bragg peaks should 
be present at n = 1/4 and n = 3/4 consistently with the 



4 



n=l/4 



n=l/2 



(a) 







□ 
















□ 








□ 
















□ 







(b) 
















c 






























X 




X 






(7i,7i) orderin: 



FIG. 3: (Color on-line) Schematic views of the plaquette (a) 
and columnar (b) phases considered in this work for electron 
densities n = 1/4, 1/2 and 3/4 (as shown). Dots, shaded pla- 
quettes and (red) thick lines correspond to electrons, singlet 
pairs resonating on a plaquette and AF bonds, respectively. 
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TABLE I: Quantum numbers of the degenerate GS of the 
various ordered phases considered here. Note that momen- 
tum (7r,0) is degenerate with (0,7r) (not shown). Standard 
notations are used to denote the irreducible representations 
of the C4v {C2v) point group (acting around point O in Fig.IIl) 
for momentum (0,0) and (tt, tt) ((tt, 0)). (*) stands for a sec- 
ond {(7r,0),^i} state. The RSPC, columnar and mixed GS 
are 2-, 4- and 8-fold degenerate respectively. 



ring exchange: 



charge ordering observed previously at the MI transition. 
Hovifever, a full characterization of the insulating phases 
is delicate since the RSPC of Fig. ^a) in fact compete 
with the columnar- like phases of Fig.[3l[b) which optimize 
the numbers of possible t2- (or dimer-) flips. For n = 1/4 
and 71 = 3/4 both RSPC and columnar phases exhibit 
(tt, 0) and (0, tt) plaquette (and charge) ordering as shown 
later. However, at n = 1/2, only the RSPC exhibits 
(tt, tt) plaquette ordering. Similarly, only the columnar 
states break rotation symmetry (for all commensurate 
densities). 

The various candidate phases hence correspond to dif- 
ferent symmetry-breaking phenomena which should be 
directly reflected in the structure of the low-energy spec- 
tra. Table |T] summarizes the quantum numbers of the 
quasi-degenerate states which should collapse onto the 
GS in the thermodynamic limit. Indeed, from linear com- 
binations of all the equivalent translated and/or rotated 
patterns obtained from each of the representative ones 
in Fig. [21 one can easily build the (orthogonal) states 
carrying the appropriate quantum numbers reported in 
Table HI Analyzing the low-energy energy spectra should 
then give clear signals on the nature of the symmetry- 
breaking phase. 



E. Analysis of the low-energy spectra: hints for 
symmetry-breaking GS 

Prior to the detailed analysis of the low-energy spec- 
tra it is useful to carry on the connection to the usual 
Quantum Dimer Model (QDM)*? and, following Ref. [H, 
add to Tt! a "diagonal" term corresponding to the 4-site 



+ Q-2Sfc-S,jnfcnKl-n,)(l-nj)] , (5) 



the summation is again over the empty squares s — 
{i,k,j,l) and the notations are the same as in Eq. ([3]). 
The Hamiltonian at the Rohksar-Kivelson point (i.e. for 
W = t2)^^ is a sum of projectors. Although our original 
problem corresponds to W = 0, analyzing the evolution 
of the system from large negative W value (for which a 
columnar phase is stabilized) up to the RK point offers 
new insights to characterize the W = phase (we shall 
not consider here the case W > t2 where localized states 
arc stabilized). This strategy has been applied success- 
fully in Ref. 12 for n — 1/2 and we extend it here to the 
more complex n — 1/4 and n — 3/4 cases. Fig. Ushows 
the low-energy spectra vs W/t2 for all three commensu- 
rate electron densities considered here. For n = 1/2 and 
n = 3/4, increasing W from negative value upwards, a 
level crossing is found around W/t2 = Wc — —1.2 sig- 
naling a (possibly first-order) phase transition. Indeed, 
our results in Fig. IHb,c) compared to the GS degener- 
acy of the ordered states listed in Table |T] strongly sug- 
gest that the columnar phase (the RSPC) is stable for 
W/t2 < Wc {W/t2 > Wc) since a group of two or four 
quasi-degenerate states (with exactly the expected quan- 
tum numbers) separated from the rest of the spectrum by 
a small gap can clearly be identified. In contrast, no level 
crossing occurs for n = 1 /4 so that no conclusion can be 
drawn yet for this case (the 8-fold degenerate "mixed" 
state of Table |T] not being completely excluded at this 
point). For n = 1/2 and n = 3/4, the "physical points" 
W = are clearly located in the region of stability of the 
RSPC. 
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FIG. 4: (Color on-line) Lowest excitations (w.r.t. the GS en- 
ergy) vs W/t2 for n = 1/4, (a) 1/2 (b) and 3/4 (c) obtained by 
ED of a 32-site cluster. Different symbols correspond to differ- 
ent irreducible representations as shown in the legend. Other 
excitations (-1- symbols in (c)) have significantly higher ener- 
gies. In (b) and (c) the vertical arrows indicate the level cross- 
ings between the {(0,0), Bi} and {{■K,n),Ai} states charac- 
teristic of the columnar state and the RSPC respectively. 
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FIG. 5: (Color on-line) Plaquette-exchange structure factor in 
the GS of the t-J-V model for n = 1/4 and J = 0.2 obtained 
by ED of a 32-site cluster, (a) Structure factor at momen- 
tum q — (7r,0) or (0, tt) vs V/\t\; (b,c) Structure factor as a 
function of momentum following the path in the BZ shown in 
Fig. [T] Open (dashed) symbols correspond to positive (neg- 
ative) t. The points labelled as "V » 1" are obtained with 
the effective model with W — 0. 



F. OfF-diagonal plaquette-exchange correlation 
function 

In order to get complementary fingerprints associated 
to a given type of symmetry-breaking and more insights 
on the nature of the insulating phases, it is of inter- 
est to compute the plaquette-exchange structure factor 
Cpi(q) = X^R^xp (jq ■ R)Cpi(|R|), the Fourier transform 
of the spatial plaquette correlations at distance R — |R|, 

CAR) = ]^E{(^nW^aG^')) (6) 

s.s' 

' (Pn(=s))(Pn(s'))}, 

where the sum is restricted to plaquettes such that 
I Us — R-s I = R and we have again subtracted the discon- 
nected part. Despite the intrinsic complexity of comput- 
ing a 8-fermion correlator, the implementation is made 
easier by the fact that the operator to average has the 
full translation and rotation symmetry of the lattice. 

We have first studied the emergence of superstructure 
peaks in the plaquette structure factor on the original t- 
J-V model itself and the results at 1 /8-filling are shown 
in Fig. [5l For V = the correlations have no structure 
in momentum space. However, for increasing V, super- 
structure peaks at momenta (7r,0) and (0,7r) (as well as 
harmonics at (tt, tt)) clearly develop. The behavior of the 
(tt, 0) peak with V follows quite faithfully the behavior of 
the "charge" Bragg peak of Fig. [2] indicating that charge 
and plaquette orderings seem to be tightly connected. 

Next, we turn to the case of the effective model describ- 
ing the insulating phase and results are shown in Fig. [S] 
At quarter-filling (n — 1/2), the large peak emerging at 




FMKrrMKrrMK r 



FIG. 6: (Color on-line) Plaquette-exchange structure factors 
in the GS of the effective t2 — J — W model vs momentum q 
along the path shown in Fig. [1] and for different choices of W 
and J as indicated in the legend (in units of t2). (a), (b) and 
(c) correspond to commensurate fillings n — 0.25, n — 0.5 
and n — 0.75, respectively. Results are obtained by ED of a 
32-site cluster. 



(tt, tt) is directly linked to the existence of the RSPC (at 
W = 0) while it is suppressed in the columnar phase at 
sufficiently negative W values, W/t2 < Wc, where transla- 
tion symmetry is restored. On the contrary, for n = 1/4 
and 3/4, since both columnar and RSPC phases break 
translation symmetry superstructure peaks at momenta 
(tt, 0), (0, tt) and (tt, tt) are always present irrespective of 
the value of W. 
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FIG. 7: (Color on-line) _Bi-plaquette structure factors in the 
GS of the effective t2 — J — W model vs momentum q along 
the path shown in Fig. [T] Same parameters and notations as 
Fig.H 

G. Diagonal Bi-plaquette correlation function: 
hints for rotational symmetry breaking 

The analysis of the low-energy spectrum strongly sug- 
gests a transition from the columnar phase to the RSPC 
by varying W, at least for n = 1/2 and n = 3/4. We 
therefore expect to be able to detect directly the signal of 
the columnar order via the emergence of some long-range 
correlations. Note that the above investigation of the 
plaquette-exchange structure factor, by detecting lattice 
translation symmetry breaking, was able to discriminate 
between columnar and RSPC phases only for n = 1/2. 
Since the rotational symmetry-breaking in the columnar 
phase is achieved via the collapse of a {q = 0, Bi} state 
onto the GS in the thermodynamic limit, it is natural to 
define a diagonal i?i-plaquette correlator as, 

s,s' 

with the diagonal plaquette operator defined as, 

' (s) = UiUj - Ukm , (8) 

and the rest of the notations is the same as in Eqs. ([3]) 
and [HI Note that there is no need here to subtract the 
disconnected part since it is vanishing for obvious sym- 
metry reasons. 

The results for the related plaquette structure factor 
(defined as before) vs momentum are shown in Fig. [7] 
and clearly reveal the appearance of a strong q = peak 
(with an additional peak at (7r,7r) for n = 1/4 and 
n = 3/4) in the columnar phase i.e. for sufficiently neg- 
ative W values, W/t2 < Wc- In contrast, for W — 0, 
which is the case relevant to the original t-J-V model, 
the -Bi-plaquette structure factor remains very flat in 
momentum space consistently with the existence of a 
rotationally-symmetric RSPC for all three values of the 
densities studied here. This clearly confirms our previous 



results for n = 1/2 and n = 3/4 and suggests strongly 
that the case n = 1/4 is similar. We therefore attribute 
the lack of level crossing in the previous spectral analysis 
for n = 1/4 to a finite size effect. 

H. Comparison to the bosonic case 

It is of interest to compare our results to the bosonic 
system (again with infinite on-site repulsion and large 
NN repulsion) whose effective Hamiltonian is given by ^ 
where -Pn(s) is now replaced by -Pn(s) = fejfejbfefo; in terms 
of the usual hard-core boson creation and annihilation 
operators on the empty square (i,j,k,l). For n = 1/2 
the bosonic system can be mapped on the usual fully- 
packed quantum loop model (or a quantized six-vertex 
model) with only kinetic processes (i.e. W = 0) exhibit- 
ing a plaquette phase^ bearing some similarities with the 
fermionic RSPG^^. However, the stability of the plaque- 
tte phase w.r.t. additional parameters (e.g. W) is more 
fragile for bosons. The n = 1/4 and n = 3/4 cases reduce, 
for bosons, to the usual QDM^ studied recently by exten- 
sive QMC methodsi^ and believed to have a columnar GS 
in the absence of diagonal interactions (i.e. W = 0). This 
is different from our results for spinful fermions exhibiting 
a rotationally-invariant RSPC GS. This emphasizes the 
special role of the spin degrees of freedom in stabilizing 
plaquette phases. Note however that Ref. 5 reports for 
the hard-core bosons system at n = 1/4, increasing the 
NN V, a MI transition at V/t ~ 3.2 towards a plaquette 
phase different from the supposed GS in the V/t » 1 
limit described by the QDM^^. The evolution from the 
MI point to larger V needs certainly to be clarified. 

I. Conclusions 

To summarize, using a t-J model extended with NN 
Coulomb repulsion, the metal-insulator transition has 
been investigated on the planar pyrochlore frustrated lat- 
tice at 1/8-filling (i.e. n = 1/4), a tractable case nu- 
merically. For increasing repulsion, the emergence of su- 
perstructure (charge) Bragg peaks signals the onset of 
the insulating phase which also immediately develops 
some form of plaquette ordering. To understand fur- 
ther the complex nature of the insulating phase, an effec- 
tive model valid for commensurate fillings like n — 1/4, 
n = 1/2 and n = 3/4 and in the limit of large NN re- 
pulsion is used. Remarkably, the frustrated nature of 
the lattice leads then to ice rule-like constraints so that 
the Hilbert space can be mapped onto a manifold of fully- 
packed dimcr or loop coverings with extra spin degrees of 
freedom. Quantum fluctuations take the form of simple 
(second-order) kinetic processes that bear strong resem- 
blance with the usual bosonic QDM. Among the various 
type of lattice symmetry-breaking states, several com- 
pelling evidences are given in favor of a plaquette GS, 
the RSPC. This new state of matter introduced earlier 
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in the quarter-filled casein, is extended here to 1/8 and 
3/8 fillings for which it has a twice larger supercell and 
a twice larger GS degeneracy. A detailed analysis of the 
low-energy spectra supplemented by symmetry consider- 
ations as well as the calculation of several types of plaque- 
tte correlations support our claims. To carry on the anal- 
ogy with the QDM, following Ref. we also introduced 
a diagonal plaquette interaction and showed signatures 
of a (presumably first order) transition towards a colum- 
nar phase as this parameter is tuned to more negative 
values. 
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